Angular dependence of novel magnetic quantum oscillations in a 
quasi-two-dimensional multiband Fermi liquid with impurities 
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The semiclassical Lifshitz-Kosevich-type description is given for the angular dependence of quantum 
oscillations with combination frequencies in a multiband quasi- two- dimensional Fermi liquid with a 
constant number of electrons. The analytical expressions are found for the Dingle, thermal, spin, 
and amplitude (Yamaji) reduction factors of the novel combination harmonics, where the latter two 
strongly oscillate with the direction of the field. At the "magic" angles those factors reduce to 
the purely two-dimensional expressions given earlier. The combination harmonics are suppressed in 
the presence of the non-quantized ("background") states, and they decay exponentially faster with 
temperature and/or disorder compared to the standard harmonics, providing an additional tool for 
electronic structure determination. The theory is applied to Sr2Ru04. 
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The low-dimensional organic conductors exhibit very 
unusual electronic properties, like magnetic field induced 
superconductivity, and are the focus of solid state re- 
search Magnetic oscillations of susceptibility and 
resistivity are the main source of detailed information 
about their electronic structure, yet there is no complete 
theory of these phenomena for quasi-two-dimensional 
(quasi-2D) metals. The analytical semiclassical (Lifshitz- 
Kosevich type Q) theory of the de Haas- van Alphen com- 
bination harmonics in the multiband two-dimensional 
(2D) metals has been suggested recently Even earlier 
it was shown Q for a system like Sr2Ru04 ||], that the 
magnetic quantum oscillations in a multiband 2D metal 
with a fixed electron density [canonical ensemble (CE)] 
are qualitatively different from those in an open system 
where the chemical potential is fixed [grand canonical 
ensemble (GCE)]. The chemical potential oscillates with 
the field in closed systems and this provides a mechanism 
for different bands to communicate with each other in 
CE producing a dHvA signal with the "sum" [Q and the 
"difference" combination frequencies [^-|| in addition to 
the ordinary dHvA frequencies. Those novel oscillations 
have been observed in quantum well structures |^ and 
Sr2Ru04 [|l^,|l^. Their amplitudes are comparable with 
the standard components, and they are robust with re- 
spect to any background (non-quantized) density of states 
at low temperatures M, but fall off exponentially faster 
with the temperature [pi. 

Many layered systems, including Sr2Ru04, are actu- 
ally quasi-2D. Weak coupling between conducting layers, 
which introduces a dispersion of the Fermi surface in the 
direction perpendicular to the layers, leads to a strong de- 
pendence of magnetic oscillations on the angle 9 between 
the normal to the conducting planes and the magnetic 
field B [n3|. This is related to the fact that the energy 



spectrum in the field becomes almost 2D (the width of 
the Landau minibands almost vanishes) at some "magic" 
(Yamaji) angles where the amplitude of the oscilla- 
tions is strongly enhanced. The theory [^ applies only to 
the clean systems at the magic angles. The full analytical 
semiclassical theory of the quasi-2D multiband canonical 
metals at finite temperature with full account for the an- 
gular dependence of the Landau miniband widths and the 
spin factors, the collision broadening of the minibands 
(Dingle factor), and the background density of states is 
given in the present paper and applied to a generic case 
of Sr2Ru04. 

First, we derive a convenient expression for a multi- 
band two-dimensional thermodynamic potential in mag- 
netic field B, mainly in the units ?i = |e| = c = 1, 
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Consider the multiband quasi-2D system with a general 
dispersion law 
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where kx,ky are the momenta in the conducting plane, 
ta is the hopping between the (conducting) layers, d the 
distance between the layers, j the integer number. We 
are mainly interested in the situation where the kinetic 
energy of electrons is much larger than the hopping be- 
tween layers, ji — ^ 2ta- The density of states (DOS) 
with an account for collision broadening of the Landau 
levels (Dingle factor) (l^Jl^ can be written as 
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where A5 = 2TT\e\B/fic is the cross-sectional area in It- 
space between the two successive Landau orbits, 

ta^nk, = + ^a{n + 1/2) + COS jk^d + gaCrHBB, 

(5) 

the energy dispersion for a present geometry of the Lan- 
dau orbits, where Da = 2taJo {jk fad tan 6) with Jq{x) 
the zeroth-order Bessel function, Ta = 7r/2TQ, Tq the 
scattering mean free time in zero field in the band 
a, and the background (bg) is included in the DOS (^) 
too 1^. After integration we obtain an important exact 
result 
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where k'j^/2ma = M ^ 

ea,n = Aqo +t^Q(n + 1/2) +g„CT^BB, (8) 

= |i? COS 0/mQ,c| the cyclotron frequency with the cy- 
clotron mass nia, A^g the band edge in zero magnetic 
field, /Lt the chemical potential, ga the electron g-factor, 
a = ±1/2, fiB the Bohr magneton. The band index 
a = ba includes the band index 6 and the spin index 
a. There pa is the zero-field density of states in the 
band a. For the energies of interest, |e — Aq| > 2ta, 
Pa = ma/2'Kh'^d per energy and unit volume, otherwise 
Pa = (rna/2iT^h^d) arccos[(e — Aa)/2ta]. Note that in a 
clean system at the magic angles, where Da = = 0, 
the expression for v reduces to i^(e — ea^n) = S {e — ea,n) , 
meaning that the spectrum becomes 2D, the Landau 
minibands reduce to the Landau levels, and the previ- 
ous expressions Q fully apply. The background density 
of states, pbg, corresponding to possible non-quantized 
(or largely broadened) bands, can be included in (||) as 
the formal limit ujbg — > 0. Obviously, the non-quantized 
background will contribute to the non-oscillating charac- 
teristics of the system. In a clean limit the density of 
states takes a standard form with the one-dimensional 
square-root singularities Af ^ l/^^D"^ — — ^a,nY [0- 
By applying the Poisson formula ||^ to the sum over n 
in the thermodynamic potential 
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with ^a = — Aq, and Aq 
written as 
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Aao + gaCr^isB, it can be 



n = no + n, (lo) 

where, after substituting x — {e — ea,n)/Da, 

dx 
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is the "classical" part of the thermodynamic potential, 
with F^ = Fa/Da- In the GCE fto does not oscillate as 
a function of 1/B, and contains the contribution due to 
spin susceptibility (Pauli paramagnetism). At low tem- 
peratures one finds 



a,bg 
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The oscillating part of the thermodynamic potential 
includes only the Landau quantized bands. 
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This expression, after integrating over e and x, reduces 
to 
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where Fb = {p ~ Abo) \mb\/ {2pBTncos9) = hcS fb/2'!Te, 
Sfb is the mean Fermi surface zero- field cross-section (in 
standard units), and the spin- related phase ipa = CTb, 
7h = 9b /{2mcos9). If the scattering time were not 
dependent on the spin projection, the phase would be the 
only quantity explicitly depending on spin in Eq. (p^), 
hence one can perform a summation over cr in the second 
term to reveal the standard spin reduction factor cos irrjb 

The amplitudes of the Fourier harmonics in Eq. ( |14| ) 
are explicitly given by 

Al = SiAur C-^) C-^^) RnR., (15, 
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where Rt (z) ^ z/ sinhz is the usual temperature reduc- 
tion factor 1^, Ry{z) — Joiz) the orientation (Yamaji) 
factor noticed in Ru = g-'^^rFb/uJb Dingle expo- 
nential damping due to collision broadening of the Lan- 
dau minibands, R^ = cos i:r^b the spin reduction factor, 
which all strongly depend on the orientation of the field. 
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FIG. 1. The amplitudes of magnetic susceptibility har- 
monics Xa [panel (c)], X/3, Xj, and X/3±c« [panel (a)] in 
Sr2Ru04 for T = Td ~ 0. Parameters used in Eq.(5) are 
d = 6.37A, j = 2 for a, j = 1 for /3 and 7 bands [10]. The 
corresponding reduction factors Ry and Ra are shown in pan- 
els (b) for P and (d) for a bands. The arrows mark the dips 
in the harmonic Xa, corresponding to zeros of the Yamaji Ry 
(Y) and the spin Rs (S) factors. The maximal amplitude of 
a, P, and /3 ± q bands is found at the first Yamaji angle 33° 
(see text). 



angle 9 (Fig. 1). As usual, one can view the Dingle factor 
as corresponding to the effective temperature T + Tu in a 
clean system, where the Dingle temperature Td = Th/n. 
For the magic angles, where Da — 0, one has Ry — 1 
and hence a pure 2D situation is recovered [^. However, 
generally in a multiband case one should expect that only 
a spectrum of a particular band will become 2D at some 
particular magnetic field tilt angle. 

The expression for tl, Eg. (p^ , contains the (first) term, 
responsible for the Landau diamagnetism and the oscil- 
lating term, which is responsible for the de Haas-van 
Alphen effect. It is small compared with the "classical" 
part, since Q/flo ^ (w/^)^ ^ 1 for the present semiclas- 
sical situation. The Fourier components appear with the 
frequencies rFa- We recover from fig (|l2|) and the first 
term in Cl (|l|) the Pauh xp = jfJ-i Y.a Pa9a and the 
Landau xl = — -jMb cos^ Pa(m/mQ,)^ susceptibili- 
ties, respectively. Since usually ga = 2, those are almost 
standard, apart from the cos^ 9 factor characteristic of 
the present geometry with the tilted magnetic field. 

Note that the chemical potential (and carrier density) 
oscillates in a closed system and, unlike in GCE, the 
"classical" part of contributes to oscillations as well 
[^ . The relevant thermodynamic potential of the closed 
system (CE) is the free energy, F = + fiN, for a fixed 
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tential is 01 



N 



-dVl/d^i. The chemical po- 



M = Mo + A, 

N _ N 

P Pq+ Pbg 



M = 



(16) 
(17) 



where /io is the non-oscillating, while /i is the oscillating 
part of the chemical potential, N = —dil/dfi is the os- 
cillating part of the density of electrons, p = J2a Pa = 
Pq + Pbg is the total density of states, including the quan- 
tized pq and the background pbg DOS. Substituting this 
expression into Qq, Eq. (p^, we obtain F = Fq + F, 
where the oscillating part is M 



F = n , 

2p' 
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while Fq is the non-oscillating part. One sees that the 
difference between the free energies of the ensembles is 
directly proportional to the fiuctuation of the particle 
density, as it should, and is suppressed by the background 
density of states pbg, ^ N'^/{pq + Pbg)- The oscillating 
part of the particle density is small, but both terms in 
( [1^ ) give comparable contribution to the magnetization 
oscillations. Indeed, N = Xia-^a; and at low tempera- 
tures 
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< 1, 



T < UJa 
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since F^ ^ B, and it falls off exponentially with tem- 
perature, Na/No ~ (r/^)exp (-27r2r/wa) , at T > coa- 
Since the density oscillations are small, sometimes they 
can be screened by the back electrode in e.g. quantum 
well structures [^ . In a more explicit form one obtains 

^ = ^ E p^-i + 4 E E cos 2.. (§ - i) 

a b r—1 ^ ■'^ 
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It is the last term, which yields the combination Fourier 
harmonics with the frequencies F = rF^ ± r'Fi,'. Their 
amplitudes. 
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are comparable with the standard single-band harmonics 
at low temperatures, T < ti)Q/27r^r, as found earlier jsj 
and confirmed experimentally pj-p^. The combination 
harmonics are suppressed in presence of the background 
density of states. The spin factor depends on the tilt 
angle, and this results in a strong angular dependence of 
the corresponding dHvA amplitudes [0,0 . Incidentally, 
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if one of the angular dependent factors vanishes for some 
particular harmonic, ro-Fb, it would notuux up with other 
bands to produce combination harmonics, whereas rFb 
(r 7^ ro) would (cf. numerical results [p^t). 

Similar to the usual 3D situation, the oscillations pro- 
duce the (partial) contributions to magnetization. Ma, 
and susceptibility, Xa, which are much larger than the 
non-oscillating contributions. Indeed, the ratio of corre- 
sponding amplitudes at low temperatures is |Xal/xo ~ 
(Fa/Bf {pam/7:^p\ma\r) Jq (2^ri?„/w„) e-2--r„/-„ » 
1, since Fa/B ^ 1. Proportionality to (Fa/B)'^ is the 
property of the two-dimensional geometry The am- 
plitudes of the standard and the combination harmonics 
can be easily found from the expressions given above. 
The ratio of susceptibilities is 



^r-F^,±r-'F^, Stt^ Tr' A^,' COS TTr' 7fc' f F 



pUJbUJb' 



(22) 



(Fb±Fb> 
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where denotes the amplitude of the standard dHvA 
harmonic for an open system, with x'^ the amplitude of 
the novel combination harmonics. We see that the com- 
bination harmonics are suppressed exponentially with re- 
spect to temperature and/or disorder compared to the 
standard harmonics. Fig. 2. In addition, they contain 
the product of the angular dependent reduction factors 
and, therefore, vary rapidly with the angle 6, their fre- 
quency being dominated by the band with the largest ra- 
tio 2tb/LiJb (e.g. band in /3±Q! harmonics in Sr2Ru04, 
Fig. 1). 

It is important to note that even in a simplest case 
of a one-band system the standard harmonics rFb are 
modified by the second quadratic term, therefore 



r-l 



p—1 p—1 
XrFb±r'F^, = 'XrFb,±r'Fy ' ^ 7^ ^'7 (25) 

and the wavcfront of magnetization of the main har- 
monics in the closed system with the carrier density 
N = const is inverted compared to the open system, 
/I — const, as it should j^,^,^. 

The present formalism is applied to Sr2Ru04 with the 
results shown in Figs. 1 and 2, calculated with the pa- 
rameters from Refs and the g-factor gb — 2 for all 
the bands. The main maxima in the Fourier harmonics of 
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with F — rFb ± r' Fb' . Thus, for r = r' = 1 we have, for 
a warped 2-band cylindrical Fermi surface without the 
background DOS, 

Xn,±F,, ^TT^rub' rJo(27r£)6Va;b')cos7r7fc'e'2'^r'''/'^'>' 9. 



FIG. 2. The amplitudes of the harmonics in Fig. 1, in- 
cluding also 2a, versus the temperature in a clean (left panel) 
and a disordered (right panel) Sr2Ru04 at the (Yamaji) angle 
e = 33°. The Dingle temperature To = 0.63K [10]. 



susceptibility x for the band a are found at 0^1(2) = 33° 
(55°), for the band (3 at ^ 9a, and for band 7 at 
71(2) — 27° (49°), which is in fair agreement with the 
data 9ai 
and 9ai = 



30.6°, 9pi = 30° and 9^i = 15.3°, Ref. 0, 



30°, 9 



/31(2) 



26° (56°), Ref. M. Those 



ima in the harmonics apparently coincide with the points 
where the Yamaji factor is unity, Ry = 1, and the system 
effectively becomes 2D, Figs. 1(b), (d). The minima, ob- 
served in a band at 24° Q (25° ^) and 40°, are very 
close to the points were both the Yamaji Ry and the spin 
Rg factors vanish: Ry = at 23° and 41°, and Rg = 
at 21° and 43°, Fig. 1(d). 

The conspicuous rapid variation with the field an- 
gle 9 of the main, XfSi Xji the combination, xp-a, 
Xp-a, harmonics. Fig. 1(a), is explained by the large 
factor 2tb/ujb in the argument of the Bessel function 
in the Yamaji factor, Eqs.(15), (p3|): 2tf}/ujp = 7.5, 
2t.j/Ljj.j = 5. This necessarily leads to the high fre- 
quency variations of the amplitudes, and apparently only 
the envelope of those oscillations has been sampled ex- 
perimentally in Refs. (see also a numerical study 
(H). Those rapid variations with angle 9 might be observ- 
able, since the condition for the Yamaji approximation, 
T:kj^{tb/eff < AS" |]l3| seems to hold for the /3-band 
[7rfcy^(f^/e/)^ : AS" « 1 : 5], while for the 7-band the 
ratio is only about 1:3. Interestingly, the spin-factor Rg 
defines the envelope of the /?— amplitude, cf. Figs. 1(a) 
and 1(b), so the studies of the minima and maxima on 
the X (9) [|lO| should allow for accurate determination of 
the g- factors gb- 

The combination harmonics contain the extra tempera- 
ture and the Dingle reduction factors, so they are falling 
off with cither temperature or disorder, or both, faster 
than the standard harmonics do, see Eqs.(21 ),(na). Fig. 2. 
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The temperature dependence of the combination har- 
monics generally cannot be characterized by some effec- 
tive mass mrb,±r'b', although at higher temperatures the 
relation mrb,±r'b' ~ rnib + r'mb> holds approximately. 
All the harmonics decay quickly with T, especially in the 
presence of even weak disorder (the Dingle temperature 
Td = 0.63K [0), Fig. 2(b), which is in very good agree- 
ment with experiment (cf. Fig. 2 in Ref. |^). 

Importantly, in the first experiments on Sr2Ru04 the 
disorder was actually larger, since the elastic mean free 
path was I ~ lO^A 0, compared to I ~2050-5000A in the 
later experiments • The higher disorder in the first 

samples, and relatively small magnitude of the combina- 
tion peaks at the field B directed along the c-axis, 6* = 0, 
instead of the "magic" angle 6 = 30° , where all the am- 
plitudes are enhanced. Fig. 1, has possibly prevented the 
discovery of the combination harmonics, predicted in Ref. 

in the first experiments on Sr2Ru04 The com- 
bination harmonics have been detected in Refs. 
The dependence of the dHvA combination amplitudes 
on angle, temperature, and disorder provide additional 
valuable tool for studying the band structure and carrier 
densities in the multiband quasi-two-dimensional metals. 
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